Effect of phase noise on useful quantum correlations in Bose Josephson junctions 
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In a two-mode Bose Josephson junction the dynamics induced by a sudden quench of the tun- 
nel amplitude leads to the periodic formation of entangled states. For instance, squeezed states 
are formed at short times and macroscopic superpositions of phase states at later times. The two 
modes of the junction can be viewed as the two arms of an interferometer; use of entangled states 
allows to perform atom interferometry beyond the classical limit. Decoherence due to the presence 
of noise degrades the quantum correlations between the atoms, thus reducing phase sensitivity of the 
interferometer. We consider the noise induced by stochastic fluctuations of the energies of the two 
modes of the junction. We analyze its effect on squeezed states and macroscopic superpositions and 
study quantitatively the amount of quantum correlations which can be used to enhance the phase 
sensitivity with respect to the classical limit. To this aim we compute the squeezing parameter 
and the quantum Fisher information during the quenched dynamics. For moderate noise intensi- 
ties we show that these useful quantum correlations increase on time scales beyond the squeezing 
regime. This suggests multicomponent superpositions as interesting candidates for high-precision 
atom interferometry. 

PACS numbers: 03.75.-b,03.75.Mn 



I. INTRODUCTION 

Confined ultracold atomic gases are promising candi- 
dates for implementing quantum simulators and for ap- 
plications in quantum technology, due to the high con- 
trollability of the experimental parameters such as the 
atomic interactions Q and the geometry of the trap . 
Among the applications we cite high-sensitivity atom in- 
terferometry, vifhich can be used for enhancing the pre- 
cision in atomic clocks and in magnetic field sensors [3- 
0]. Of particular interest are Bose Josephson junctions 
formed by two modes of a Bose-Einstein condensate. The 
modes may correspond either to two internal states of 
the condensed atoms in a single potential well or to two 
spatially separated wavefunctions in a double well. Dur- 
ing the dynamics following a sudden quench of the tunnel 
amplitude connecting the two modes, squeezed states are 
formed at early times. It has been shown theoretically 
[Toj and experimentally that these states can be used 
to estimate phase shifts with sensitivity below the shot 
noise limit, the limit one obtains using classical states. 
The highest possible phase sensitivity, limited by quan- 
tum uncertainty only, can be achieved by using macro- 
scopic superpositions of e.g. atomic phase states [Tll.[l^. 
Such superpositions arc however formed at later times 
during the quenched dynamics of the BJJ [l3l - [T5| . They 
are expected to be very fragile with respect to decoher- 



ence effects caused by particle losses [16[ , collisions with 
thermal atoms [ill j interaction with the electromag- 
netic field Il9i . and random fluctuations of the trapping 
potential [20|. 



In this work we consider the effects of phase noise on 
the states formed during the quenched dynamics of the 
BJJ. Phase noise is induced by stochastic fluctuations of 
the energies of the two modes of the BJJ. Together with 
atom losses, such a noise is one of the main sources of 
decoherence in the experiments of Refs. 0, 0, HH. In- 
terestingly, macroscopic superpositions of phase states 
in BJJs have been shown to be robust with respect to 
phase noise, their decoherence rate being independent of 
the total number of atoms in the condensate [l^] . Be- 
low we show that these long-lived states can be useful 
in interferometry to improve phase sensitivity. In par- 
ticular, we compare the best possible phase sensitivity 
obtained with the state of the BJJ at the times of forma- 
tion of macroscopic superpositions to the one obtained at 
earlier times when squeezed states are produced. This al- 
lows us to determine which are the most useful quantum 
states for interferometric applications in the presence of 
phase noise. In order to quantify the amount of quan- 
tum correlations useful for interferometry, we calculate 
the quantum Fisher information. In the theory of esti- 
mation of an unknown parameter, this quantity is related 
to the bound on the precision with which the unknown 
parameter - in interferometry, the phase shift - can be 
determined [13, [l^ ■ We show that for moderate noise 
intensities the quantum Fisher information at the time of 
formation of the first superpositions of phase states ex- 
ceeds the one found at the time at which squeezed states 
appear. In other words, despite the action of phase noise, 
a better phase sensitivity can be reached by using states 
produced at times later than the time for optimal squeez- 
ing. 
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The paper is organized as follows. In Sec. II we recall 
the definition of the two parameters relevant in interfer- 
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FIG. 1. (Color online) Rotations on the Bloch sphere in the 
interferometric scheme: the input coherent state at the north 
pole (green disk) is rotated around the y-axis by an angle 7r/2 
(blue disk) and afterwards around the z-axis by the unkown 
phase ip (black disk). The precision Atp on the estimation 
of ip is larger than the size \/N /2 of the disk, representing 
the angular momentum fluctuations, divided by the radius 
N/2 of the sphere. The last rotation around the y-sxis is not 
represented. 

ometry, i.e. the coherent spin squeezing and the quan- 
tum Fisher information, and link them to multiparticle 
entanglement. In Sec. Ill we present the model which de- 
scribes the quenched time evolution of a BJJ both in the 
absence and in the presence of noise. A peculiarity of the 
quenched time evolution is the formation of multicompo- 
nent superpositions of phase states, and we illustrate the 
effect of phase noise on those states in Sec. IV by calcu- 
lating a suitable probability distribution. In Sec.V we 
compute the coherent spin squeezing and the quantum 
Fisher information during the quenched dynamics of the 
BJJ, first in the absence and then in the presence of noise, 
and study quantitatively the loss of useful correlations as 
a function of time. Finally, Sec. VI contains some con- 
cluding remarks. 



II. PHASE ESTIMATION IN ATOM 
INTERFEROMETRY 

A. Phase estimation and the quantum Fisher 
information 

The goal in interfcromctry is to estimate an unknown 
phase shift (p with the highest possible precision. In atom 
interfcromctry, an input state is first transformed into a 
superposition of two modes, analogous to the two arms of 
an optical interferometer. These modes acquire distinct 
phases Lpi and (^2 during the subsequent quantum evolu- 
tion. They are are finally recombined to read out interfer- 
ence fringes, from which the phase difference = v^i — (^2 
is inferred. The interferometric sequence can be de- 
scribed by means of rotation matrices acting on the two- 
mode vector state, that is, by SU(2) rotation matrices in 



the Schwinger representation [23 . [25| . The generators of 
the rotations are the angular- momentum operators Jx, 
Jy, and Jz, related to the annihilation operator aj of an 
atom in the mode j = 1,2 by Jx = (a|a2 + a2ai)/2, 
Jy = —i{a\a2 — 0201)72, and Jz = h = {a\ai — 0,^0,2) /2, 
the latter being the number imbalance operator. Let us 
consider the case where the two modes correspond to two 
internal states of the atoms in an optically trapped Bose- 
Einstein condensate. The total number N of atoms in the 
condensate is assumed to be fixed and all atoms are ini- 
tially in the lower energy state (mode j = 1). The input 
state is then \nz = N/2), where jn^ = n) = \ni,n2) de- 
notes the Fock state satisfying Jz\nz = n) = n\nz ~ n), 
Til = N/2 + n and n2 = N/2 — n being the number of 
atoms in the lower and upper modes, respectively. The 
application of a 7r/2 laser pulse with frequency in res- 
onance with the two internal levels plays the role of a 
beam splitter in optical interferometers. It brings the in- 
put state onto the coherent state \d = tt/2^ 4) = 0), where 
the SU(2) coherent states are defined as [2g 

with Q = tan(6'/2) cxp(— i^i). It is easy to show that 
|6l,0) cx (e-"^sin(6l/2)4 -I- cos(6'/2)4)^|0) (here |0) is 
the vacuum state), meaning that all atoms occupy the 
same one-atom state. A coherent state can be visualized 
as a disc of diameter \/N /2 on the Bloch sphere of radius 
-/V/2, centered at iV(sin6' cos (p, sin 6* sine/), — cos 0)/2. The 
coordinates of the center are the expectation values of the 
angular momentum operators Jx, Jy, and Jz m |0, 0), 
whereas the diameter of the disc gives the quantum fiuc- 
tuations of Js = J • n in the directions n tangential to the 
sphere. The coherent states with 6* = 7r/2 on the equator 
of the Bloch sphere are referred to as phase states. The 
Fock state jn^ = A^/2) is a coherent state with 6 = tt, lo- 
cated at the north pole of the Bloch sphere. The action 
of the beam splitter is therefore a rotation of the atomic 
state around the y-axis by an angle of n/2 radians, lead- 
ing to the phase state \9 = 7r/2,(/) = 0). Then the state 
is rotated around the z-axis by the free evolution, the 
phase accumulation being due to a different energy shift 
between the two states. This rotation is the analog of 
the different phase paths in the two arms of an optical 
interferometer. The consecutive rotations of the input 
state on the Bloch sphere are represented in Fig[TJ Fi- 
nally, by recombining the two paths, the state is rotated 
again around the y-axis by an angle of —Tr/2 radians. 
The interferometric sequence can thus be described by a 
succession of three rotations, and the output state of the 
linear interferometer is 

IV^out) = e-'i^ye-''^J^e'iJy\iji^) = e-'^^-'^l^i^) , (2) 

where IV'in) is the input state, assumed here to be pure. 
More generally, the output state of the interferometer is 

Pontic) =e-''^J^Pine'^J^\ (3) 
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where pin is the input density matrix and n the unit 
vector representing the effective rotation axis associated 
to a given interferometric sequence. 

In a typical experiment one has access to the probabil- 
ity distribution associated to the operator Jz measured 
with respect to the output state. This quantum distri- 
bution depends on the phase shift ip. The latter is then 
determined by means of a statistical estimator depending 
on the results of the measurements of Jz in the output 
state. The precision Atp with which the phase shift (p 
can be determined depends on the chosen estimator, on 
the input state and on the measurement performed on the 
output state. Optimizing over all possible measurements, 
the best precision that can be achieved for a given input 
state Pin is, according to the Cramer- Rao bound [23| . 



Aip > (A.^)best 



1 









Pint 'Jn 



(4) 



where m is the number of measurements and Fglpin, Jn] 
the quantum Fisher information given by (23j 



Fc 



Pin; 'Jfi 



= 2 E 



jPl -Pm)'^ 
Pi + Pm 



\JnH\\ (5) 



{|/)} being an orthonormal basis diagonalizing pin = 
(with pi > Q and Y.1P1 = !)• The Fisher 
information ([5|) depends on the input state and on the 
direction fi of the interferometer. It has the meaning 
of the square of a " statistical speed" at which the state 
evolves along the curve defined by Eq. ^ in the smce of 
density matrices when the parameter ip is varied [ill [23j : 
if one increases p starting from p = Q with a fixed ve- 
locity the faster the state ([3]) becomes distinguishable 
from Pin, the larger is its quantum Fisher information 
Fq . Hence the bound ([5]) relates the problem of estimat- 
ing a phase shift in an interferometer to the problem of 
distinguishing neighbouring quantum states |23|. Indeed, 
the quantum Fisher information is related to the Bures 
riemannian distance on the space of density matrices (27j . 

For pure input states IV'in), the quantum Fisher infor- 
mation is given by the quantum fluctuation (AJ,-j)^ = 

(■0in|i||V'in) - (V'inlinlV'in)^ of Js, 



Fq \i^in)Jn =4(AJ3) 



(6) 



This allows to reinterpret the Cramer-Rao lower bound 
(|4]) as a generalized uncertainty principle 



A(/? A J,-j > 



1 



(7) 



in which the generator j^j of the transformation ([3]) and 
the phase shift ip play the role of two conjugate vari- 
ables - (p being here not an observable but a parameter 
[23} . For instance, the Fisher information of a phase 
state \'4'in) = \9 = 7r/2, 0) in the directions n = ex. 



By, and 6*2 are equal to iVsin^0, N cos^ (j), and N, re- 
spectively. According to ([7]), for such a state the best 
precision that can be achieved on the phase shift is 
(Aiy9)bcst = 1/vNm = {Aip)sN, corresponding to the 
shot-noise limit of independent atoms. 

The saturation of the bound ^ requires both a suit- 
able classical post-processing on the m outcomes of the 
measurements (e.g. the maximum likelihood estimation 
in the limit of large m and the knowledge of the 

optimum observable to measure. This latter task can be 
difficult as the optimum measurement may depend on 
the phase shift itself psl ]. 

It can be shown [ill , [l^ that for any separable input 

state Pin, Fglpin, Jn] < N, SO that 



Fq Pin, Jfi > N 



(8) 



is a sufficient condition for pin to be entangled; in other 
words, Fq — is an entanglement witness. By Eq.(|4]), 
the inequality ([5]) is a necessary and sufficient condition 
for sub-shot noise sensitivity (A(/?)bcst < iA(p)sN- In 
what follows, the input states leading to such a condi- 
tion are called useful states for interferometry (or, more 
briefly, "useful states"). It is worthwhile to stress that 
the inequality ([8]) is not a necessary condition for en- 
tanglement: indeed, there exists entangled states which 
are not useful for interferometry, that is, with a Fisher 
information Fq < iV [III [H ■ 

The quantum Fisher information is bounded by iV^. 
This is easy to show for pure states by noticing that the 
largest square fluctuation of J^ in Eq.® is smaller or 
equal to (see [l3l); for mixed states this follows 

from the convexity of Fq (see [ll|). According to Eq.(|4]), 
the best sensitivity that can be achieved in linear inter- 
ferometers [131 is then (A(^)bost = {Alp)hl = 1/A^. This 
corresponds to the so-called Heisenherg limit. This limit 
is reached using highly entangled atoms as input state, 
e.g. the macroscopic superposition given by the so-called 
NOON state Ii/'noon) = (|A^, 0) + e*"|0, N))/V2, with a 
being a real phase. The quantum Fisher information of 
a NOON state is equal to A^^ in the direction n = Cz- It 
is instructive to compare this result with the value of the 
quantum Fisher information for a statistical mixture of 
the same states, pnono = {\N ,Q){N ,Q\ + \Q,N){Q, N\) /2. 
The latter is found with the help of Eq.® to be equal 
to A^ in all directions n in the (xOj/)-plane and to van- 
ish in the direction e*j. Therefore, the scaling of Fq 
like N'^ for pnoon = |'0noon)('0noon| is due to the 
presence of the off-diagonal terms /Onoon — Pnono = 
{e-'°'\N, 0)(0, A^l + e*"|0, N){N, 0|)/2. 

To summarize, the study of the quantum Fisher infor- 
mation and its scaling with the number of atoms allows 
to quantify the amount of quantum correlations which 
can be used to enhance the precision on the phase shift 
in interferometry. 
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Phase estimation 


Entanglement 


Fq [pin] > iV (Av?)bcst < (A(^)SJV 


Fq [Pin] >N^pinJ^ Pscp 


[pin] < 1 ^ (A(^)bost < {A(p)sN 


[/5in] < 1 => Pin / pscp 



TABLE I. Necessary and/or sufficient conditions for sub-shot 
noise phase sensitivity in an atom interferometer and multi- 
particle entanglement in terms of the quantum Fisher infor- 
mation and spin-squeezing parameter. 



B. Coherent spin squeezing 

Atomic squeezed states are examples of nonclassical 
states useful for interferometry, which have been recently 
realized experimentally 043 ■ The coherent spin squeez- 
ing parameter quantifies the aiigular- momentum fluctu- 
ations along the direction n d, according to 



Sn 



Pin-) Jfi 



(9) 



where 



pi = cos (pCx + sm By 

P2 = — cos 6 sin (pCx + cos 9 cos e^, + sin 6 (10) 



are unit vectors perpendicular to 



n = sm a sm oCx ~ sm f cos m Cy 



cos a ez, 



(11) 



and (•) = tr(-pin) is the mean expectation in state pin. A 
state Pin is said to be squeezed in the direction n if the 
squeezing parameter satisfies 



P'uii Jfi 



< 1. 



(12) 



It is known that Eq. (fT2|) provides both a sufficient (but 
not necessary) condition for sub-shot noise sensitivity [lo| 
and a sufficient (but not necessary) condition for entan- 
glement of Pin ■ We remark that the squeezing criterion 
([T2|) does not recognize all useful states in interferome- 
try. For instance, the NOON state does not fulfill this 
criterium even though it leads to the best achievable pre- 
cision. The criteria for entanglement and sub-shot noise 
sensitivity are summarized in Table L 



C. Optimum coherent spin squeezing and quantum 
Fisher information 



The quantum Fisher information .Fq and the spin 
squeezing parameter ^ introduced in the previous sub- 
sections depend on the direction of the generator which 
defines the intcrfcrometric sequence ([3]). For instance, as 
shown in Sec lII Al Fq[\iPnoon) , Jz] = N^, corresponding 
to a maximally entangled state, whereas in the perpen- 
dicular directions FqI^jj^oon), Jx] = ^q[|^noon), -'^y] = 
N. Hence, in order to quantify the useful correlations of 



a quantum state, one needs to optimize Fq and ^ over 
all the possible directions by defining [g^ 



[pin] = min ^1 [pin] , Fq [pin] = max Fq 



(13) 

Let us consider the 3x3 real symmetric covariance matrix 
7[Pin] with matrix elements 



lij [Pin] 



E 



jPl - Prnf 
Pi +Pm 



SRe 



{l\J^\m){m\Jj\l) 
(14) 

where {|/)} is the orthonormal eigenbasis of pin as in 
Eq.([5]). According to standard linear algebra, the maxi- 
mum of FQ[pin, Jfi] — 4(n, 7 [pin] n) over all unit vectors 
n is equal to 



Fq [pin] = 47n 



(15) 



7max being the largest eigenvalue of the matrix 7[/5in]- In 
the following it will be useful to define also the matrix 

G,, [p] ^ i {J J, + J, J,) - (J,) (J,), (16) 

where (...) ~ tr(. . . /5), with p being the system density 
matrix. Note that for pure input states \ipin) the matrix 
7ij[|'(/'in)] reduces to the matrix G^j [|V'in) (■'/'in|]i which is 
easier to compute than the more general expression (|14p . 
The optimum quantum Fisher information is then given 
(up to a factor four) by the largest uncertainty of the an- 
gular momentum operators (see Eq.([6])). For the sake 
of brevity, in the following we will omit both the adjec- 
tive "optimum" and the explicit dependence on the input 
state, designating the optimum coherent spin squeezing 
and the optimum quantum Fisher information respec- 
tively by and -Fq, unless where source of confusion. 



III. QUENCHED DYNAMICS OF A BJJ 

A. Noiseless dynamics 

We describe a Bose Josephson junction (BJJ) by a two- 
mode Hamiltonian (Slj . which in terms of the angular- 
momentum operators introduced in SecHIl reads 



H 



(0) _ 



XJ' - XJz - 2KJx 



(17) 



This Hamiltonian models both a single-component Bose 
gas trapped in a double- well potential [2l| - external BJJ 
- and a binary mixture of atoms in distinct hyperfine 
states trapped in a single well [l,!!^! - internal BJJ. In the 
external BJJ the two modes i correspond to the lowest- 
energy spatial modes in each well. For the internal BJJ, 
the two relevant modes are the two hyperfine states. The 
first term in the Hamiltonian p7p describes the repulsive 
atom-atom interactions; for the external BJJ, x is half 
the sum of the interaction energies Ui in the two modes. 
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whereas for the internal BJJ x = {Ui + f^2)/2 — C/12 also 
depends on the inter-species interaction C/12. In both 
cases, A = + {N ~ l)(C/2 ~Ui)/2 is related to the 
difference /S.E = E2 — Ei between the energies of the two 
modes. The last term in (jl7p corresponds to tminelling 
between the two wells or, for the internal BJJ, to a reso- 
nant laser field coupling the two hyperfine states, which 
can serve to implement a 50% beam splitter as described 
in Sec lII Al Both x and K are experimentally tunable 
parameters. In internal BJJs arbitrary rotations of the 
form ([3]) can be performed with a suitable combination 
of laser pulses. Such rotations are typically realized fast 
enough to neglect the non-linear effects induced by the 
interactions Q . The residual effect of interactions on the 
interferometric sequence has been recently addressed in 
Refs. [ll,!!!. 

We consider the dynamical evolution induced by a sud- 
den quench of the tunnel amplitude K to zero, taking as 
initial state the phase state \9 = 7r/2,(^ — 0). This is 
the ground state of the Hamiltonian (fT7|) in the regime 
KN ^ X where tunelling dominates interactions (in the 
internal BJJ this state can be produced by using a laser 
pulse as explained in Sec lII Ap . We assume a fixed total 
number of atoms N, that is, we do not account for atom 
losses. Going to the rotating frame [ssf . we may suppose 
that A = 0. In the absence of noise, the atomic state 



V2,0-O) 



(18) 



displays a periodic evolution with period T = 27r/x if 
N is even and T/2 if N is odd, corresponding to the 
revival time. At intermediate times the dynamics drives 
the system first into a squeezed state at short times, then 
at times tq = 7r/(xg) to a macroscopic superposition of q 



phase states given by |13l-ll 



9-1 
UO ^ Ck,q 
k=0 



,<Pk,q 



(19) 



where |mo| - q-^'\ <l>k.q = (2fc - N)^, Ck,q = e-'^"/" 
if q is even, and (pk^q = (2fc + 1 - iV)^, c^,, = 
^ink{k+i)/q if g ig rpj^ig follows from Eqs.dl]) and 

([T8|) and the use of the Fourier expansion e~"^"i/^ = 
'^o ELo e*'^'^'/«e-2*'^"i'=/«. As in the case of the NOON 
state discussed in Sec lII A[ the two-component superpo- 
sition formed at t = i2 = T/4 leads to the best achiev- 
able phase sensitivity if used as an input state of the 
interferometer described by Eq.Q [lH. We will show 
in SecHI] that the multicomponent superpositions with 
3 < g ^ "v/iV, which are formed at earlier times tq < t2, 
lead to comparable phase sensitivities up to a factor of 
two. The fact that squeezed states and macroscopic su- 
perpositions of phase states are intrinsically produced 
by interatomic interactions yields a major advantage of 
atomic interferometers over optical ones. 



B. Dynamics of a noisy BJJ 

The presence of noise during the dynamical evolution 
of the BJJ affects the premration of the aforementioned 
useful entangled states [331 • We focus here on phase 
noise caused by a randomly fluctuating energy difference 
AE{t) between the two modes, assuming that the inter- 
action energies Ui and U2 are not fluctuating. In the 
single- well experiment (internal BJJ), such a noise is 
induced by fluctuations of the magnetic field, whereas in 
the double-well experiment [2l| (external BJJ) it is in- 
duced by fluctuations in direction of the laser beam pro- 
ducing the double- well potential with respect to the trap- 
ping potential. The corresponding time evolution which 
follows the sudden quench iiT — > is described by the 
time-dependent Hamiltonian 



H{t) = xl -X{t)J, 



(20) 



Even in presence of noise, the time evolution following 
the quench can be exactly integrated since the noise term 

MtjJz commutes with the noiseless Hamiltonian x^z 
[2^ . For a given realization of the stochastic process 
A(t), the state of the atoms at time t is 



|^(f)) = e-"^<^'^J'\iP<^°\t)) 



(21) 



where (j){t) = — drXlr) and IV'^-^HO) is the time- 
evolved state (|18p in the absence of noise. The system 
density matrix is then obtained by p{t) = \ip{t)){'ip(t)\ = 
/ dP[X] \ip{t)){il;{t)\, where the overline denotes the av- 
erage over the noise realizations. The introduction of the 
distribution probability for the random angle 0(t), 



/(0,t) = J dP[A(t)]<5(0-0(t)) (22) 
allows to write it as 

d</)/(0,t)e-'^^"^p(o)(t)e'^^"% (23) 



where p'--^\t) = l^p'^^^t)) (ip'^^^t)] is the density matrix in 
the absence of noise. Under the hypothesis of a gaussian 
noise (see AppendijjX| the probability distribution ([^^ 
reads 



2TTa{t) 



(■» + At)^ 



(24) 



where A = AE + (N - l)(C/2 - Ui)/2 and the variance 
{t) is given in terms of the noise correlation function 



h{T-T') = X{t)X{t')^X = AE{t)AE{t')~AE (25) 
by 



a^{t) ^ dr dT'h{T~T'). 
Jo Jo 



(26) 
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Wc note that h depends only on the time difference t — t' 
by the stationarity of the stochastic process A(t), which 
also implies X{t) = A(0) = A; moreover, h decreases to 
zero at sufficiently long times. By projecting Eg. ([^5]) on 
the Fock basis {jn^ = n)} we obtain 

K = n\p{t)\n, = n') = e" 

x(n,=n|p(")(t)|n, =n'). (27) 

In order to discuss the effect of the phase noise on the 
state of the atoms we briefly discuss the noise variance 
a{t). We first notice that under our hypothesis, a{t) and 
thus the decoherence factor (given by the first exponen- 
tial in the right-hand side of Eq.([27])) is independent of 
the number of atoms N in the BJJ. This is in contrast 
with the usual scenario for decoherence which predicts 
stronger decoherence as the number of particles in the 
system is increased. As a consequence of this fact, macro- 
scopic superpositions of the form (fT9|l are robust against 
phase noise, as was shown in [2^ and will be detailed in 
SecHV] below. 

Let us denote by tc the largest time such that /i(t) ~ 
h{0) = dX{0)^ = (5A^ and by Tc the characteristic time 
at which /i(r) vanishes. If the time evolution occurs on 
a short scale such that t < tc then the colored nature of 
the noise plays an important role and 

a^{t) criSX^t"^. (28) 

If instead the time evolution occurs on a time scale much 
larger than the noise correlation time Tc we obtain the 
same result as for white noise, 

/•OO 

a^{t) ~ 2t / h{y)dy, (29) 
Jo 

which corresponds to the Markov approximation. 

The effect of phase noise can be part ially suppressed 
by using the so-called spin-echo protocol JSg . This strat- 
egy was followed in a recent experiment [6| . The analysis 
discussed in this section can be adapted to take into ac- 
count the residual effect of phase noise when spin echo 
pulses are applied, see Appendix iBl 



A. Structure of the density matrix in the Fock 
basis 

In the absence of noise the quenched dynamics 
of the Bose Joscphson junction leads to the forma- 
tion of coherent superpositions with q components as 
given by Eq. p9p . The corresponding density matrix 
p(o)(t^) = |V;W(i,))(V'(°Hi9)l has the form p(°\t,) = 
Sfe k' Pk'k'i^i)^ where the indices k and k' label the var- 
ious components of the superposition and p^^i},{tq) = 
q~^Ck,qcl,jTT/2,(j)k,q){n/2,(j)k',q\. For general decoher- 
ence processes one expects that, by increasing the in- 
tensity of the noise, p'-^'>{tq) will evolve into the statis- 
tical mixture of phase states Pkki^i)^ moreover, the 
larger the atom number N the weaker should be the noise 
strength at which this occurs [s^, H^l ■ It was found in 
[2^ that for the phase noise considered in Sec lIIII the ac- 
tual scenario for decoherence is different from the usual 
one. Indeed, the typical noise intensity at which the co- 
herences between distinct phase states |7r/2, (t)k,q) are lost 
turns out to be independent of the atom number. This is 
a consequence of the fact that the decoherence factor a{t) 
is independent of A^, as shown in Sec lIII Bl Furthermore, 
for superpositions with a large number of components 
this intensity is larger than the noise intensity at which 
phase relaxation occurs. In what follows we discuss the 
origin of this fact. 

Since the noise is expected to destroy correlations be- 
tween different components, we decompose the density 
matrix in its diagonal (intra-component) and off-diagonal 
(intercomponent) parts, focussing on the latter one to 
quantify the decoherence. We have then p^°^ = p^^^ + 
where 

p'-^\t,) = 'EPkkit,) (30) 

k=0 

and 

P'flit,)= E P^'(t,)- (31) 

k,k'=0;k^k' 



IV. EFFECT OF PHASE NOISE ON 
MULTICOMPONENT MACROSCOPIC 
SUPERPOSITIONS OF PHASE STATES 

Before analyzing the Fisher information and the spin 
squeezing during the quenched dynamics discussed above 
in detail, we wish to study the nature of the state of the 
atoms under phase noise at the specific times tq which in 
the noiseless BJJ correspond to the formation of multi- 
component superpositions of phase states. We first illus- 
trate the effect of the noise on the structure of the density 
matrix, then we study a suitable probability distribution 
which is particularly sensible to decoherence. 



Using Eq.O and the identity e2*'=("'-")'^/9 = q ii 

n = n' modulo q and otherwise, the matrix elements of 
p''p(tq) in the Fock basis are 

{n, ^ n\p^^\tq)\n, = n') (32) 

I if n' 7^ n mod q 

where p is an integer and / denotes the integer part. By 
using p^°]{tq) = e-^^^ll'i\d = 7r/2,,/> = O)(0 = tt/2,4> = 
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0|e' 



Pd itq)i we also get 
s(0) 



{nz = n\pmtg)\n^ = n') 



(33) 



TV 



2JV V^+n; \^+n' 



N 



li n' = n + pq 
ii n' ^ n mod q 



witha„„' = {n' + n- N){n' -71)71 /q. The use of Eg. ([771) 
allows to obtain the corresponding expressions in the 
presence of noise, 



(f^|pd,od(^q)|?^') = e 



"|pSod(i.)l"'> (34) 



up to a phase factor irrelevant for decoherence, with aq = 
a{tq). In the strong noise limit aq ^ 1, the off-diagonal 
part Pod of the atom density matrix vanishes whereas the 
diagonal part pd tends to a matrix which is diagonal in 
the Fock basis, 

Pd{tq)^Poo = j ^\0 = 7r/2,^){e = 7r/2,^\ (35) 

n=-N/2 ^2 / 



The fact that the diagonal part of the atom density ma- 
trix decays faster than the off-diagonal part for increasing 
noise strengths is readily explained by examining the 
structure of the noiseless density matrices in Eqs. (|32|) and 
(pSj) . The first off-diagonal elements of Pd{tq) in the Fock 
basis are those for which n' = n zt q while the first off- 
diagonal elements of Pod(tq) satisfy n' = n± 1. Hence, it 



results from Eq. (p4| that the off-diagonal elements of pd 
vanish at the noise scale a ~ 1/q while the off-diagonal 
elements of pod vanish at the larger noise scale a ~ 1. 
In other words, the noise is more effective in letting pd 
converge to pao than in suppressing pod , and this effect is 
more pronounced the higher is the number of components 
in the superposition. An illustration of such anomalous 
decoherence is given in Fig. 2. The middle panels show 
that for intermediate noise strengths, pd has already ac- 
quired its asymptotic diagonal form psp . while pod has 
not yet vanished. As wc will see in Sec lVBI below, these 
results imply that, for moderate strengths of phase noise, 
macroscopic superpositions are formed and provide quan- 
tum correlations useful for interferometry. 




FIG. 2. (Color online) Matrix elements of the diagonal (intra- 
component) part pd{to,) (panels a),c),e)) and the off-diagonal 
(inter-component) part Pod(i3) (panels b),d),f)) of the den- 
sity matrix in the Fock basis at time t = tz as the noise is 
increased from as = (a),b)) to as = 0.9 (c),d)) and as = 2.9 
(e),f)). 



these distributions, which would be absent if the atoms 
would be in a statistical mixture of phase states. 

The probability distribution of in the state p can 
be calculated by a straightforward generalization of the 
calculation in [41| , as the Fourier coefficient of the char- 
acteristic function h^{r]) = tr[e~"''^*p], namely, 



(36) 



For the quenched dynamics of the Bose Josephson junc- 
tion in the presence of noise, the characteristic function 
reads 



B. Angular momentum distributions 

The anomalous decoherence of the atomic state can be 
visualized by plotting the probability distribution P^{r) 
of the eigenvalues of the angular momentum operators 
Jit> = Jx sin (j) — Jy cos (j) in an arbitrary direction of the 
equatorial plane of the Bloch sphere [4l| . The presence of 
correlations among the components of the superposition 
formed at time tq is indicated by interference fringes in 



N/2 

h^{r];t) = ^ gnn'{t){nz = n\p'^^\t)\nz = n) 

xDn'n{-(l>,vA) (37) 

where 5„„'(i) = ^-a\t){n-n'f/2^{\ti,i-n') g^^^ 

-Dri'n(~0, 0) is the matrix element of the rota- 
tion operator e""'"^* in the Fock basis, which is given by 
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FIG. 3. (Color online) Probability distribution P^/2(7', ^s) of 
the eigenvalues of Jx for the three-component coherent super- 
position (solid lines) at increasing noise strength from as = 
(a),b)), to as = 0.9 (c),d)) and as = 2.9 (e),f)) with = 
20 atoms. The blue dashed curves indicate the large-noise 
intensity and large A*' limit given by Ea. (|39p . 



file displays two peaks corresponding to the projections 
on the X-axis of the phase states \0 = 7r/2,(/) = (/ffc^a), 
0fc,3 = 0, 27r/3, 47r/3 (the "phase content" of the state, 
accounted for by Pd{ts,)) and interference fringes, due to 
the coherences between these phase states (contained in 
Podit?!.))- In the presence of noise (b)-c)), the phase profile 
of each component of the superposition spreads and the 
characteristic peaks in the profile of the distribution are 
smeared out (phase relaxation). At strong noise intensi- 
ties, Pd{tq) approaches the steady-state given by the den- 
sity matrix (|35p . which is symmetric in the (a;0?/)-plane. 
As a consequence, the corresponding probability distri- 
bution P^{r, oo) = P(r, oo) = tilpaoWx = r){nx = r\] 
is independent on (j). In the semi-classical limit ^ 1, 
this distribution can be easily calculated since takes 
the values cos 0/2 in the phase states \TT/2,(f>) apart 
from small relative fluctuations of the order of 1/\/N. 
Hence, recalling that poo is a statistical mixture of the 
states |7r/2,0) with equal probabilities (see Eq. ipS]) ). 



P{r, oo) = c / ( 
Jo 



N \ 1 

dcjjS [ — cos (j) — r ] ~ —- 



1 



(39) 

where c is a normalization factor. The semi-circle law 
([39ll is indicated by the blue dashed curve in panel c) of 
Fig El For finite A^, one finds 



P(r, oo) 



N 



ir[f 



r [A^ + 1] 



On the other hand, the vanishing of Pod(tq) tends to 
diminish the contrast of the fringes in the distribution 
Prf,{r,tq), until they are completely washed out in the 
asymptotic distribution (panel c) of FiglS]). The fact that 
phase relaxation occurs at a lower noise strength than 
decoherence is evident in the panel b), where the profile 
of P0(r, tq) is close to the asymptotic distribution P{r, oo) 
while interference fringes due to Pod{tq) are still visible. 



(see e.g. [H, Eq. (D6)) 



min{Ar/2-ri',Ar/2+n} 

E (-1) 

fc— max{0,n — n'} 

m 

X 



n\e "''""In^ 



-ir;J^ I 

N 



N 



f +n' 



k)\{f+n-k)lk\{k + n'-n)l 



(38) 



2k-\-n' —n 



1] 

cos — 
2 



N+n-n'-2k 



— im\n —n) 



The probability distribution in the absence of noise de- 
rived in [ilj is recovered by setting gnn'{t) = 1 in Eg. ([57)) . 

The distribution P.,/2{rM) = = r\i'^^\hW of 
the eigenvalues of Jx (satisfying Jx\nx r) ~ r\nx = r)) 
is shown in FigjS] for the three-component superposition 
of phase states in the absence of noise (panel a)). Its pro- 



V. QUANTUM FISHER INFORMATION AND 
COHERENT SPIN SQUEEZING DURING THE 
QUENCHED DYNAMICS OF A BJJ 

We present in this section the estimate of the use- 
ful quantum correlations which are formed during the 
quenched dynamics of the Bose Josephson junction intro- 
duced in Secs lIIII and IIVI For this purpose, we evaluate 
the quantum Fisher information and the coherent spin 
squeezing parameter. We consider first the noiseless evo- 
lution for which analytical expressions can be obtained, 
then we present the numerical results for the dynamical 
evolution in the presence of noise. 



A. Dynamics in the absence of noise 

In the absence of noise the atoms are in a pure state 
\^^'^\t)) during all the dynamical evolution. The co- 
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variance matrix 7*'°''(i) associated to this state is thus 
given by Eq. p6p . For the quenched dynamics described 
in SecHniby using Eqs. (P), one finds that 



(J, 



\(0) 



yit 



(40) 



(J. 



v(0) 



(t), (41) 



where (..)^'''' = tr(.../5(°) (i)) and z/'^^jt) corresponds to 
the visibility of the Ramsey fringes The angular- 

momenta covariance matrix (jl6p finally reads 







[(7V-1) 



cos 



N-2 



(2t) - (TV + 1)] 



\ m^cos^-^r)sm{T) 
where we have introduced the rcscaled time t ~ 2TTt/T = xt and 



7i°)(r)^((AJ.)^)(") 



N r 



— ^i ^ cos 



(r) sin (t) 



4 



(0). ^ N 

'± V ; 16 



(42) 



(iV - 1) cos^-2 (27-) + (AT + 1) _ 2Af cos2(^-i) (r) 

o L 

The two other eigenvalues of the matrix (|42p are 

-(Af - 1) cos^-2 (2r) + (iV + 3) ± (Af - 1)y/(cos^-2 (2r) - 1)2 + 16cos2(^-2) (^) gj^^ (r) 



(43) 



(44) 



(see also Ref.|43|). We remark that the matrix (|42|) has 
the property that its eigenvalues at times r and tt — r 
(and, similarly, at 27r — r) coincide, hence it suffices to 
discuss its behaviour at times t belonging to the interval 
[0,r/4] (i.e., re [0,7r/2]). 

According to Eq. (|15p , the quantum Fisher information 
is given by the largest eigenvalue, 



i^Q(r) =4max{7i°)(T),7f(r)} 



(45) 



We demonstrate in Appendix [C] that the coherent spin 
squeezing ^ is always optimum along a direction con- 
tained in the (?/Oz)-plane. The optimal spin squeezing 
parameter ([T3|) is thus related to the lowest eigenvalue 

7i^''(r) of the submatrix 7^'^-' (r) obtained by removing 
the first lign and column in the matrix (j42p . Using 
Eqs.lgni) and (gT]), one gets 



7Vi.(0)'(T) 



(46) 



The direction of optimum squeezing is given by the eigen- 
vector 



^(0). > ^(0)/ N 

(r) = n_ [T] 



- sin ' (r) Ey + cos of^ (t) (47) 



I 

associated to the eigenvalue ^''^^ (r) of the covariant ma- 
trix. One easily finds 



9f\r) 



= arctan 



v(0) 



Cretan ( (^'^^"^->)- 

2 \ I ?2\(0) / f2\(0) 



(48) 



where {•, •} denotes the anticommutator. 

The direction of optimization n'p of the quantum 
Fisher information is either given by (if "ix^ > 7+ "* ) or 
by the eigenvector fi^^ associated to the eigenvalue 7^"' 



(if Jx^^ < ) • The latter condition is satisfied at times 
shorter than t* , see Appendix |Dl As both these eigen- 
vectors are orthogonal to n^'' (since the matrix 7^°) is 
symmetric), it follows that coherent spin squeezing and 
quantum Fisher information are optimized in perpendic- 
ular directions. At short times, when the state of the 
system is a squeezed state, this has a clear physical inter- 
pretation: the quantum Fisher information is maximum 
in the direction of maximal angular momentum fluctu- 
ations, which is perpendicular to the direction of lowest 
fluctuations yielding the best squeezing. 

The Fisher information (|45|) and the squeezing param- 
eter (|46p obtained from Eqs. (03]) and are shown as 
a function of time in Fig. |4l At short times, the co- 
herent spin squeezing (a)) is below one, indicating the 
presence of a squeezed state. To compare spin squeez- 
ing and Fisher information, we introduce the parameter 
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FIG. 4. (Color online) a) Coherent spin squeezing and b) 
quantum Fisher information during the quenched dynamics 
of a BJJ with = 100 atoms as a function of time (in units 
of the revival time T) in the absence of noise. The dashed 
line in the second panel represents the parameter = N/^^ . 
Horizontal and vertical gridlines in panel a): minimum of 
the coherent spin squeezing and corresponding time tmin (see 
text), c) Non-optimized quantum Fisher information along 
the a;-axis (dashed line) and the y-axis (dotted line) . For com- 
parison, the optimum quantum Fisher information of panel b) 
is also shown (gray solid line). The vertical gridlines corre- 
spond from right to left to the time t = t/s of formation of 
the first macroscopic superposition, see Ea. (|49|) : to t = t* , 
see Appendix [Pl and to t = T/4 — t/s. The horizontal grid- 
line shows the shot-noise level Fq — N. d) Angles O'^^ in 

Eq. (|48I) (dashed line) and 9p^ (solid line) giving the optimiz- 
ing direction for the spin squeezing and the quantum Fisher 
information as a function of time. 



= N/^"^. This parameter was shown in [TT| to coincide 
at short times with the Fisher information, indicating 
that Fq and ^ provide essentially the same information 
for squeezed states at such times. This property is il- 
lustrated in Figl4]-b), and is demonstrated in the large 
N hmit in Appendix |El At larger times, the squeezing 
parameter first reaches a minimum and then grows to val- 
ues larger than one (that is, F^ decreases and becomes 
smaller than N). This does not imply that the atomic 
state is not useful for interferometry since, as described 
in Sec nil the squeezing criterion is only a sufficient con- 
dition for useful entanglement [ll| . Indeed, the quantum 
Fisher information increases above the shot noise level 
Fq = N until it reaches a plateau value Fq ~ A^(A^+l)/2 
(see Appendix [E|). at the time 



(49) 



tf,^T/VN 



This time corresponds to the time of formation of the 
"first" (in chronological order) multicomponcnt superpo- 
sition, as one can infer from the following argument. The 
largest number of phase states of size -\/]V /2 (see Sec lII Ap 
which can be put on the equator of the Bloch sphere of 
radius N/2 is g^ax ^ 2i:N/\/N = 2tt\/N. The time of 
formation of the multicomponcnt superposition with the 
highest number of phase states is tj^ = T/(2ginax): lead- 
ing to Eq. (|49ll . Note that tfs is also the time scale for 
phase diffusion, that is, for the decay of the visibility (|4ip . 
It is seen in FigUthat Fq diplays a sharp maximum at 
t = t2 = T/A, in correspondence to the two-component 
macroscopic superposition which has the highest possible 
Fisher information Fq = N'^. This result is not surpris- 
ing since this two-component superposition (fT9|) is the 
analog in the phase variable of a NOON state, from which 
it can be obtained by a 7r/2 rotation around the y-axis. 
Panel c) of Figl4] shows the Fisher information in the di- 
rections Cx and Cy. In the time regime corresponding 
to the plateau, they are almost equal due to the sym- 
metry of the multicomponcnt superpositions (this means 

that the eigenvalues ji^'^ and 7^^'' are almost degenerate). 
As one approaches the two-component superposition, the 
optimizing direction changes to the a;-axis, which is the 
direction of maximal angular momentum fluctuations. 

In panel d) of FigHl the angle o'^^ giving the direc- 
tion of highest spin squeezing in the {yOz) plane is rep- 
resented as a function of time together with the corre- 
sponding angle for the Fisher information, which 

gives the optimizing direction of the Fisher infor- 
mation according to Eq. pTjl . Table II summarizes the 
aforementioned results. Some analytical results obtained 
for short, intermediate and long times in the limit ^ 1 
of a large number of atoms are given in AppendislEl 
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Time 


Optimum quantum Fisher information Fq 


Optimizing direction 


t = 


N 


degenerate in (yOz) plane 


< t < T/N 


/• , given by Ea.(|ET|l 


- cos (t) ey - sin df^ (t) 


T/N < t < tn,in 


^ , given by Eg.dEl]) 


— 


^min ^ ^ ^ tfs 


^ , 3i/3jv5/3 < i?Q < 0.4323iV2 


- Sy 




— ^ , Fq ~ iV(iV + l)/2 


- Sy 


t* <t< T/4 


, N[N + 1)/2<Fq<Fq{T/4) = N^ 


Ex if N is even, ~ if A'^ is odd 


Time 


Optimum coherent spin squeezing parameter = N/^^ 


Optimizing direction 


t = 


N 


degenerate in (yOz) plane 


< t < T/N 


/- . F^^Fq 


- sin 6»f ' (i) Cy + cos e»f ' (t) e. 


T/N < t < tmin 


/• , < < F^(Wn) = 2iV5/33-2/3 


~ Cz 


^min ^ i ^ tfs 


\ , < F^ < 2iV^/^3~^''^ 




tfs«t< r 


\ , < Fe < A/- 




t* <t< T/4 


\ , < Fj < iV/3^/2-1 





TABLE II. Optimum coherent spin squeezing parameter, op- 
timum quantum Fisher information and corresponding opti- 
mizing directions during the quenched dynamics of a Bose 
Josephson junction in the absence of noise for A'^ 3> 1. The 
arrows indicate whether the function is increasing or decreas- 
ing with time in a given time interval. 



B. Dynamics in the presence of noise 

Let us now consider the effect of the phase noise intro- 
duced in Sec lIII Bj on the results obtained in the previous 
subsection. For coherent spin squeezing the calculation 
can be carried out analytically. We start with the obser- 
vation that even in the presence of noise {Jy)t ~ {Jz)t ~ 
and more generally the angular-momentum eovarianee 
matrix G defined in Ea. (|16p has the same structure as 
the matrix (|42|l in the noiseless case. Therefore, the ar- 
guments used in Appendix [C] can be taken over to the 
noiseful case. We thus conclude that the squeezing pa- 
rameter is minimum in the (?/Oz)-plane, and is given 
by Eq. (|46| , evaluated for the corresponding quantities in 
the presence of noise. In particular, the bare visibility 
i>^^\ Ea. (|4T|l . should be replaced by the visibility ly in 
the presence of noise [12] , 



N 



(J. 



(50) 



and 7^°^ by the lowest eigenvalue G_ of the restriction of 
the eovarianee matrix G to the (?/Oz)-plane. 

We are now going to compute G_ and the spin squeez- 



ing parameter explicitly. In order to do so, we need to 
perform the averages in the presence of noise using the 
full density matrix p{t): (...)* = tr{. . . p{t)). They are 
related to those in the absence of noise according to 



dcl>f{c^Me"^''U-''^'')f^ (51) 



where the expectation value inside the integral is taken 
for the pure state \ij}^^\t)) in the absence of noise. The 
rotated angular momentum operators in the above ex- 
pectation value are equal to cos (jiJx — sin (jjjy , sin (jjj^ + 
cos^Jj,, and Jz for i = x^y, and z, respectively. A simi- 
lar derivation holds for ({J,;, Jj})* = ir[p{t){Ji, Jj}], with 
the result 



I _ g-2a"(i) 



f2\(0) 



l + e-2°'(*) 



y 2 
({J„Jj)t-e-"^(*)/2({J„JJ)(°) 

{{JxJy})t^ {{JxJz})t^Q. (52) 

Finally, the submatrix matrix G'{t) reads 



G'(r) 



N ^_g-2a^(r) _ ^) ^^^^-2 (3^) + {N + 1)] 



-N{N ~ l)cos'''~^ (r) sin(T) 



■N{N~ 1) 



cos 



N-2 



(t) sin (r) 



N_ 

4 



(53) 
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Thus, by Eqs.(|46l), and one has 

1 



e-'''(^)(7V - 1) cos^-2 (2t) + e°'(^)(7V + 3) 



-{N - l)e"'(^y(l - e-2'^'(^) cos^-2 (2t))2 + 16e-'^'(^) cos2(^-2) (r) sin^ (r) 



(54) 



The angle which identifies the optimal squeezing direc- 
tion is given by Eq. (H51) , in which the matrix 7^"^^ should 
be replaced by G' . 

We proceed by illustrating our results for the squeezing 
parameter in the presence of phase noise. For the calcula- 
tions we have chosen a noise range of direct experimental 
relevance, as extracted from the fit of the visibility decay 
data in Fig. 4. 15 of Ref. with our prediction given by 
Eg. ([50]) . For the noise variance a2(r) we take the short- 
time behavior a'^{T) = {6\/xYt'^ expressed by Eq. ((28l) 
since the experimental visibility exhibits a gaussian de- 
cay even for small interactions x [35| . This indicates that 
in the time regime < t ^ t/s, the phase noise in the 
experiments has strong time correlations (colored noise). 
The squeezing parameter as a function of time is shown 
in Figl5]-a). As seen in the figure, the presence of noise 
degrades the squeezing, as its minimum value increases 
at increasing noise strength. We also notice that the time 
for optimal squeezing t^^^ is slighly shorter than in the 
noiseless case. As is shown in Appendix |E] in the limit of 
a large number of atoms we find that the minimum value 

of ^2(^) is 



Sn 



(55) 



this minimum being reached to leading order at the same 
31/67V-2/3 noiseless case. 



time Tir 



-(0) 



This means that by increasing the number of atoms the 
noise becomes less efficient in limiting the highest squeez- 
ing which can be reached during the dynamical evolution. 
This results from the fact that the time fmin at which 
this highest squeezing is produced is shorter for larger 
N , whereas the effect of the noise on the density ma- 
trix ([?7|) is independent of N , as stressed in Scc lIIIBl 
The angle 9(^{t) which identifies the optimizing squeezing 
direction is represented in dashed lines at various noise 
levels in Figl5]-d) . As discussed in Appendix [El similar 
to the noiseless case, 0^ almost vanishes at intermediate 
times 7V"i < r < 7r/2 - TV^^. 

The evaluation of the optimum quantum Fisher infor- 
mation (|14|) requires a numerical diagonalization of the 
density matrix p{t) given by Ea. (l27p . For the time de- 
pendence of a?{t) we take again the short-time approxi- 
mation given in Ea.([28ll. even if there is no experimental 
evidence which justifies such a choice at times t ^ T. 
This choice corresponds to the worst possible scenario for 
decoherence, as in the markovian regime the dependence 
of a2(f) is weaker (see Eq.i^) [H]. The behaviour of Fq 
as a function of time in the presence of noise results from 
the competition of two phenomena: (i) in the absence 



of noise, at short times the quantum Fisher information 
grows from its initial value Fq = N to the plateau value 
Fq = N{N + l)/2 in a time interval tf, - T/y/N which 
shrinks as N becomes larger, and (ii) the decoherence 
exponent (t) is independent of TV and also grows with 
time. As a result, Fq reaches a local maximum at a time 
^max ~ tsf, with a value which increases with N and 
decreases with the noise fluctation SX^. 

The quantum Fisher information as a function of time 
at various noise levels is shown in Fig. O The short- 
time evolution is similar to the one found for the noise- 
less case, the accumulation of noise correlations being 
not yet effective. In particular, one observes that Fq 
coincides with the squeezing parameter F^ = A^/^^ at 
sufficiently small times (panel c) . For not too large noise 
intensities, Fq displays a plateau at those times which 
in the noiseless BJJ correspond to the of formation of 
macroscopic superpositions. The value on the plateau is 
smaller than in the absence of noise but it is still much 
above the shot noise level Fq = N. This indicates the 
presence of useful correlations which remain in spite of 
the decoherence effects induced by the noise. This effect 
is due to the robustness of the multicomponent super- 
positions [13] with respect to phase noise discussed in 
Sec. IIVI above. For higher noise levels, the width of the 
plateau is reduced and the peak at t2 = T/i correspond- 
ing to the two-component superposition in the absence 
of noise disappears completely, meaning that decoher- 
ence has washed out the useful quantum correlations at 
t2 (three bottom curves in the Figl5]-b)). In the fimit 
of very large noise intensities the Fisher information at 
times tq of formation of g-component superpositions in 
the noiseless BJJ is degenerate in the [xOy) plane and 
tends to the asympototic value 



Fq[Po 



N{N -I) 
2N + 2 



(56) 



which can be readily obtained from Eas. (|T5)) and ((35l) . As 
illustrated in FiglH apart from short times and around 
the peak at t2, the optimization direction is in the (xOy)- 
plane and -Fq is almost degenerate in all directions of this 
plane, as in the noiseless case. 

As a partial summary, the analysis of the time evo- 
lution of the quantum Fisher information indicates the 
build-up of useful quantum correlations at times beyond 
the spin-squeezing regime. In the following we quantify 
this effect by studying the dependence of Fq with the 
noise strength and the particle number. 

Figure [7] shows FQ{t) on a logarithmic scale, evalu- 
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FIG. 5. (Color online) Coherent spin squeezing and quan- 
tum Fisher information in the presence of noise as a func- 
tion of time in units of T during the quenched dynamics of a 
BJJ. The parameters used are = 100, x = ""Hz. a) Spin 
squeezing for (from top to bottom) 5\ — 15, 10, 5, and 
OHz. Horizontal and vertical gridlines: minimum of ^'^ and 
corresponding time imin- b) Fisher information Fq for (from 
top to bottom) 5\ = 0, 0.4, 1, 2, 5, 10, and 15 Hz; the horizon- 
tal and vertical gridlines correspond to Fq = (A*' -f- 1) /2 
and t = tfs — T/\/N. c) Zoom on the quantum Fisher 
information (solid lines) and F^ = N/^"^ (dashed lines) for 
5\ = 0,2,5, 10, and 15 Hz (from top to bottom), d) Angles 
Of and giving the optimizing direction of Fq (solid lines) 
and (dashed lines) as a function t/T, for the same noise 
levels. 
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FIG. 6. (Color online) Direction-dependent quantum Fisher 
information in the presence of noise as a function of time 
in units of T during the quenched dynamics of a BJJ with 
A' = 100 atoms and x = ttHz for: a) 5\ = 2Hz, b) 
5Hz, c) lOHz and d) 15Hz, calculated along the {Ox) di- 
rection (dashed lines), the (Oy) direction (dotted lines) and 
the optimizing direction (light-gray solid line). After a time 
t ~ T/^/N (left vertical gridlines) the three values are almost 
the same, showing that the Fisher information is almost de- 
generate in the ixOy) plane, except around t = r/4 if Fq has 
a peak at this value (panel a)). The vertical and horizontal 
gridlines represent the times t — tfs and t = r/4 — tfs and 
the value of the Fisher information in the limit of large noise 
intensities given by Eq. (|56|) . 




FIG. 7. (Color online) Values of the Fisher information at 
its local maximum at time tmax (solid line, circle markers), at 
time t2 (dot-dashed line, star markers) and at the time fmin of 
maximal squeezing (long-dashed line, blue cross markers) in 
a logarithmic scale, as a function of the energy fluctuation S\ 
(in Hz). For comparison we also plot the squeezing parameter 
F^ — N/^"^ at the time tmin (dashed line, green cross mark- 
ers) in a logarithmic scale. Gridlines, from top to bottom: 
Heisenberg limit A''^ (solid), approximate value (2/3^''^) A'^/^ 
oi F^itmin) in the absence of noise, see See l V Al (dashed) . shot 
noise limit (solid), and limit of Fq for large noise intensities 
(solid) given by Ea. (|56|) . The parameters used are A' = 400 
and X = Hz. 



14 



ated at the time t = t2 = T/A of formation of the two- 
component superposition in the noiseless BJJ, as well 
as the maximum (i^Q )inax of FQ{t) in the time interval 
< t < T/8. This maximum corresponds roughly to 
the value at the plateau in Figl5l that is, to the value 
of Fq (t) at the times of formation of the first multicom- 
poncnt superpositions. It can be seen that in the range 
of noise considered (-F'Q)max stays above the shot noise 
level, and is also larger than the value FQ{tynin) at the 
time ^min of highest squeezing. The two-component su- 
perposition, formed much after the superpositions with 
a large number of components, appears to be too much 
degraded by noise to lead to any advantage in interfer- 
ometry with respect to separable states. Hence, in this 
regime multicomponent macroscopic superpositions pro- 
vide a convenient alternative to both the squeezed states 
and the two-component macroscopic superposition. 

We next study the scaling of the quantum Fisher in- 
formation with the particle number, taken at the time 
^max as before. As it is illustrated in FiglH at such a 
time Fq displays a power-law behaviour Fq ^ with 
an exponent f3 depending on the noise strength. This 
exponent is extracted from a log-linear fit of the numer- 
ical data, varying N between 50 and 400 [4J] , the latter 
value being realistic in the experiments 0. We notice 
that in the noise range considered (3 is larger or equal to 
5/3, which is the exponent corresponding to the squeezed 
state at t = imin in the absence of noise (see Sec IV Ap (isj . 
This confirms the potential improvement in interferome- 
try given by the state at time tmax with respect to the use 
of squeezeed states in the presence of phase noise. For 
comparison, we also show the scaling of Fq at the time 
t2- At that time, /3 decays faster with the noise strength, 
reaching rapidly the shot noise limit /3 = 1. This is due 
to the fact that the noise exponent a^{t) increases with 
time. 
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VI. CONCLUSIONS 

In this work we have studied the effect of phase noise 
on the formation of nonclassical states useful for interfer- 
ometry created during the quenched dynamics of a Bose 
Josephson junction. The knowledge of an exact solution 
for the dynamical evolution of the state in the presence 
of phase noise has allowed us to calculate the quantum 
Fisher information as a function of time and its scaling 
with the particle number at various noise strengths. Due 
to the anomalously slow decoherence induced by phase 
noise on macroscopic superpositions of phase states, for 
a realistic choice of noise strengths we have found that 
multicomponent superpositions are more useful for inter- 
ferometry than squeezed states. Such superpositions are 
built during the dynamical evolution of a noiseless junc- 
tion at larger times than squeezed states. The time of 
formation of the superposition with '-^ \/N components 
depends inversely on the square-root of the total number 
of atoms N . When phase noise is affecting the unitary 



FIG. 8. (Color online) a) Quantum Fisher information eval- 
uated at the time of its local maximum tmax (blue solid line) 
and at the time t2 (blue dashed line) as a function of the 
number of particles A'^ for 5A = 15 Hz, as compared to the 
shot noise limit (black solid line). Panels b),c),d),e): same as 
in a) in a semi-logarithmic scale, for various noise strengths 
5\ — 2,5, 10, and 15 Hz (from left to right and top to bot- 
tom), f) Exponent /3, extracted by a log-linear fit of the data 
in a), as a function of the energy fluctuations SX (in Hz) for 
t = imax (solid line, circle markers) and for t = t2 (dot-dashed 
line, star markers). We used x = ""Hz. 



dynamics of the junction, these multicomponent super- 
positions therefore provide an interesting alternative to 
the use of the more popular two-component superposi- 
tion, which would appear later in a noiseless junction. 



15 



ACKNOWLEDGMENTS 

We thank C. Gross and A. Smcrzi for useful discus- 
sions, and P. Hyllus for suggesting us Ref-H^]- We ac- 
knowledge financial support from CNRS, the MIDAS 
project, the Handy-Q project, and the project ANR-09- 
BLAN-0098-01. 



Appendix A: Derivation of the probability 
distribution t) for gaussian noise 

The probability distribution /((/), t) of the angle (j){t) = 
— /q dTX{T) is defined as an average over the noise real- 
izations induced by the functional P[\{t)] as in Eq. (P^ . 
or, by Fourier expansion, 

f{(t),t) = — f dP[X{t)] f due-^-^^'^e*"^". (Al) 
27r J J_oo 

We a re left to evaluate the Fourier transform of the av- 
erage e-™'^(*) = / dP[X{t)] e-™'^^*). This is readily done 
under the hypotesis of a gaussian noise. 



g-V(0(t)-0(t))2g-m</,(t) 
-■^ /o drfg dT'h{T-T') -iuM 



g-^a (t)^-iu\t 



(A2) 



where /i(r — r') = A(r)A(r') — A(0) is the noise corre- 
lation function and we used that for gaussian variables 

with ^ = one has e™? 



Using Ea. (|A2|) wc 
obtain the expression (j24p in the main text, according to 



27r 



(A3) 



Appendix B: Partial suppression of phase noise by 
spin-echo pulses 

In a recent experiment 0, phase noise was partially 
suppressed by a spin-echo protocol [s^. Let us assume 
that the state of interest (for instance, a squeezed state 
in [^) is produced after an evolution time t under the 
Hamiltonian ([20]) . In the spin-echo protocol, two short tt- 
pulses are send by a laser in resonance with the energies of 
the two modes at times t/2 and t. The effect of these laser 
pulses is to reverse the direction of Jz , mapped into — Jz , 
in the evolution between t/2 and t. Since the noiseless 
part of the Hamiltonian ([20]) is quadratic in J^, it is not 
affected by the pulses, while the noise part is linear in Jz 
and is reversed after half of the evolution. This allows to 
suppress the effect of the noise if it is strongly correlated 
between the two time intervals [0,i/2] and [i/2, tjjwhich 
appears to be the case in the experiment of Rcf . Q (see 
also dH). 



Our model in Scc lIIIB] can be easily adapted to take 
into account the residual effect of phase noise when the 
spin-echo pulses are applied. The derivation follows the 
same lines as in the main text. Ea. (P7)) still holds pro- 
vided that we use (/)(t) = /g dr Sgn(T — t/2)A(T), with the 
sign function defined as Sgn(2;) = ±1 for ±x- > 0. This 
leads to 

(Bl) 

We focus on the short time regime t < tc- The approx- 
imation /i(t) ~ h{0) yields no contribution to a^^i^oi^)- 
An expansion to second order is needed, /i(r) = h{0) + 
Th'{Q) + T^h"{0)/2 + 0(t3). Using the parity /i(-r) = 
/i(t) of the correlation function (which implies h'{0) = 0), 
we obtain [ivf 



oit) 
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(B2) 



Comparing Eqs. (|B2p and ([28|) . one sees that the effect 
of the noise at times t < is considerably reduced with 
respect to the case in absence of spin echo. 



Appendix C: Demonstration of Ea. (|46p for the spin 
squeezing parameter 

In the following we show that the spin squeezing pa- 
rameter £,'^{t) in a Bose Josephson junction is always op- 
timized along a direction contained in the (y02;)-plane. 

Let us observe that the angular momentum covariance 
matrix G{t) defined by Eq. (jl6p has vanishing matrix el- 
ements Gxy{t) = Gxz{t) = 0. In fact, in the absence 



of noise this matrix G(t) 



'{t) is given by Ea. (|42)) . 
and we have seen in Sec. IV Bl that it preserves the same 
structure in the presence of a phase noise. Thanks to this 
special structure of G'(i), the fluctuations of the angular 
momentum operator along an arbitrary direction n given 
by Eq.([ll|) is 

AJfl(t) - n,G,j{t)nj (CI) 

= s\T? 6 sir? tpGxxit) + niGij{t)nj . 



The sum over i,j in the second line can be written as 
(sin^ 6 cos^ (p+cos'^ 6)n'^ G' {t)n' , where we introduced the 
notation G'(t) for the two-by-two submatrix of G(t) in 
the plane (yOz) and the normalized vector 



y/sh?~9 



COS^ 



cos^ 9 



(C2) 



in this plane. Furthermore, we observe that during the 
dynamics of the noisy junction one has {Jy)t = {Jz)t = 
at all times. As a consequence, the expectation values 
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of the angular momentum operators along the directions 
defined by Eq. (fTO|) are given by 



{Jpi)t = COS(j){Jx)t 

{Jp2)t = ~ COS sin (j){X)t . 



(C3) 



Combining these results and using the faet that Gxx (t) > 
0, we obtain from Eq.® 



sin^ Osin^ 4>Gxx{t) 
1 — sin^ (J) sin^ 9 



> G-{t) = ^ niin {n'^G'{t)n} 

n' , \\fi' II —1 

{ffG{t)n] 



(C4) 



mm 

ne(yOz),\\n\\ = l 



where v{t) = 2{Jx)t/N is the visibility and G_(i) the 
smallest eigenvalue of G'{t). Since it is clear that the in- 
equality in Eq. (jC4p is an equality for n equal to the corre- 
sponding eigenvector n-{t) of G-{t), this demonstrates 
that the squeezing is minimized along a direction n-{t) 
contained in the (2/Oz)-plane. Combining Eas. (jl3p and 
(jC4p . we obtain that the optimum coherent spin squeez- 
ing is given by Ea. (|15)) . 



Appendix D: Determination of the time t* when the 
optimization direction of the Fisher information 
changes in the absence of noise 

If the number of atoms is even, the direction of opti- 
mization n'p} of the Fisher information in a noiseless Bose 
Josephson junction is along x-axis at the time ^2 = T /A 
of formation of the superposition of the two phase states 
\e = 7r/2,0 = 0) and \e = 7r/2,'/> = tt). These phase 
states are indeed diametrically opposite on the equator of 



Bloch sphere along this axis. Since n^^(r) = itj^' {t) is in 
the (j/Oz)-plane at times r = 2-Kt/T < 1 (see Sec lVAp . 
the optimizing direction thus changes abruptly from the 
(yOz)-plane to the x-axis at some time r* <s]0, 7r/2[ sat- 
isfying 



?(0)/ 



(0)/ *N 

7+ (t" ) 



In this appendix we determine r* explicitely in the limit 
of large total atom number A^, supposed to be even. We 
may infer from the previous discussion that r* is nei- 
ther close to nor close to 7r/2. Consequently, wc look 
for a solution of the implicit equation (jPip in the in- 
terval T G [iV~",7r/2 — A^~"], a being a positive ex- 
ponent strictly smaller than 1/2. Introducing the vari- 
ables u = cos(t) g [0, cos(A^~")] and v = cos(2t) G 
[-cos(2iV-"),cos(2iV-")], we obtain with the help of 
Eqs.d^ and (gj) 



4(7f (r)-7r(r)) 



(0)/ 



N 



-{N -ly^^^ + Nu 



,2N-2 



(D2) 



+2(iV - l)u'''-\l - u') + 0(iVu*^^-«) + 0{Nv'''-^) 



Setting ^^\t) = jx'' (j) gives the equation 

N-2 



JO), 



2-\ 



2-^^^^+0{e-^'-'^). (D3) 



For large TV, the right-hand side of Ea. (|D3p is strictly 
larger than one and is of the order of unity. Hence the 



solution must satisfy |2 



> 1 and 2 



~ ±1. 

We may exclude the positive sign as the values u = ±1 
correspond to r ~ or r = tt outside the studied time 
interval. The relevant solution u of Ea. (|D3p is thus close 
to l/\/3 and smaller than this number. Let us note that 
for odd iV's, such a solution does not exist; indeed, in 
this case Ea. ()Dip has no solution (see Sec IV Ap . Let us 
set u = l/(\/3(l + 5)). Then from Eq.(|D3]) we obtain 



g(JV-2)I„(l+65+0(5^)) ^ 5 ^ o(<S) + O (^1^ 



from which we find 



In terms of the dimensionless time r* we get 
1 \ ln(5/3) „ f 1 



V3/ 6V2A^ 



O 



(D4) 



(D5) 



(D6) 



Appendix E: Analytical results for the 
spin-squeezing parameter and the quantum Fisher 
information in the large-Af limit 

Short time regime. At times shorter than the time 
of formation of the first macroscopic superpositions, ie 
< t < t/. (or T < 1/V^), one has 7^"^(i) > -i'i\t). 
A short-time expansion in Ea. (|44|) followed by the large 
N limit yields 



F^\T)=A^X>{r)c.N 



(0), 



2^2 



N't 



-Nt\ I 



(El) 

In this time regime, the visibility (|4ip is almost equal 
(Dl) to one. In order to compare Fq with ~ A^/?^, we 



determine the ratio 
F^'\r) 



N + 1 N-1 



Ff(r) 



2 



■ cos 



N-2 



(2r) 



(E2) 



-(N-lfcos'^^-^T) sin2(r) 

by employing Eq.(|46p and identifying the product 

7+ "* (t)!-'' (t) with the determinant of the 2x2 subma- 

trix 7(°)'(i) in Eq.^. At short times t < N^^/^, which 
a posteriori turns out to be the time of optimal squeez- 
ing, the RHS of Ea.(|E2p can be approximated by one, 
yielding 

^ ^^"\r) (E3) 



17 



as discussed in the main text. At later times t N 
the right-hand side of E c|.(| E2l) can be approximated by 
1 + iV^rVe. We obtain 



whereas the squeezing parameter 

Ff^(T)~7Vz.(")'(T) 



(E7) 



N{l + N'^T^/6) 



(E4) 



where Fq\t) is given by Eq. (|Eip up to a relative correc- 
tion 0{N~^/^). The squeezing parameter (jE4[) reaches 
a minimum (C^in)^ - {3/N)^^^/2 at the time rj^l = 
3i/6jY-2/3 limit ^ 1, as assumed above (see 

also Ref. [1] where a different definition of ^ is used, which 
however almost coincides with ours at times t <^ tsf be- 
cause ~ 1). The value of the Fisher informa- 
tion at r = T^^l is FQ(r^°|j ~ The direc- 
tion of optimization of Fq is in the (yOz)-plane and is 

given by the eigenvector rt^''(r) orthogonal to 7i^''(r), 



that is, cI>P{t) = and = ^^^'^^ + ""Z^- 

finds by using the first equality in P5|) that tan 6'^'''' (r) ~ 
{Nt/2 + ^i + Ar2^2/4)-i_ The angle 61^°^ starts from 
7r/4 at T = and quickly decreases to the value 0, to 
which it is almost equal at times r ^ N~^. At such 
times ^ and Fq are optimal along and ey, respectively 
(see Eq. (|47|) ). These results are summarized in Table II. 

Intermediate times. In the time regime 5t < t < T/A — 
St with St 3> tfs, the covariance matrix (j42p takes the 
simple following form in the limit A^ ^ 1 







^N{N + l) 



(E5) 



decreases with time as i/'-^^t) decreases (second panel 
of FigU]). We have shown in Appendix [Cl that if A^ is 

even, the optimizing direction fi^-'(r) of the Fisher in- 
formation changes as r increases from the (yOz)-plane 
to the X-axis at the time r* ~ arccos(l/-\/3) defined by 
"fx'\T*) = 7^''(t*). Note, however, that any direction in 
the (a;Oy)-plane gives a Fisher information almost equal 
to the optimized value A^(A^ -I- 1)/2, as mentioned above 
and as it is clear from the structure of the matrix (jESp . 
For an odd number of atoms N, the optimal direction 
n^'' (r) remains in the (2/Oz)-plane all the way up to 



T = 7r/2 (i.e., 7i°V) < ^r any r G [0,7r/2]). 

More precisely, it is almost along the y-axis (which is the 
symmetry axis of the superposition (jl9[) formed at t = 
at times A^^^ < r < 7r/2. 

Times t close to t2 ~ T/A. At times t such that \t2 — 
t\ <^ T, the Fisher information is given by 



Ml 



N 



A^ + l + (A^-l)e-2^('^/2-r)= 



(E8) 



It increases monotonously from the plateau value (|E6p at 
times t c:^ T/A — tsf to the value A^^ at the time t = t2 oi 
formation of the two-component macroscopic superposi- 
tion, which has the highest Fisher information Fq ~ N"^ 
allowed by the Heisenberg bound. The optimal direction 
of Fq is along the x-axis if A^ is even and the y-axis if A^ 
is odd, and that of ^ is along the z-axis in both cases. 

Results in the presence of noise. The formula general- 
izing Ea. (jE4p for small nonzero noise intensities a(T) ^ 
A^"^ reads 



1 + A^Vye + NiSX/xfr"^ + o{N-^/^) 



(E9) 



Hence the Fisher information has a plateau at the value 



(E6) 



where we assumed r < A^^^/a ^^^^^^ SX/x N^^^. The 
minimum value of ^(t)^ is given by Ea. (|55p of the main 
text. The angle 9^ which identifies the optimal squeezing 

direction satisfies tan((?j) = e""^'*^/^ tan(6'^°^). 
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